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ABSTRACT 

Assuming that giant planets are formed in thin protoplanetary discs, a "3-D" system - i.e. a 
planetary system composed of two (or more) planets, whose orbital planes have large values 
of mutual inclination - can form, provided that the mutual inclination is excited by some 
dynamical mechanism. Resonant interactions and close planetary encounters are thought to 
be the primary inclination-excitation mechanisms, resulting in a resonant or non-resonant 
system, respectively. If by the end of planet formation the system is dynamically "hot", then a 
phase of planet-planet scattering can be expected; however, this need not be the case in every 
system. Here we propose an alternative formation scenario, starting from a system composed 
of three giant planets in a nearly co-planar configuration. As was recently shown for the case 
of the solar system, planetary migration in the gas disc (Type II migration) can force the 
planets to become trapped in a multiply-resonant state (similar to the Laplace resonance in 
the Galilean satellites). We simulate this process, assuming different values for the planetary 
masses and mass ratios. We show that, such a triple resonance generally becomes unstable, as 
the resonance excites the eccentricities of all planets, and planet-planet scattering sets in. One 
of the three planets is typically ejected from the system, leaving behind a dynamically "hot" 
(but stable) two-planets configuration. The resulting two-planet systems typically have large 
values of semi-major axes ratio (a - 01/02 < 0.3), while the mutual inclination can be as high 
as 70°, with a median of ~ 30°. These values are quite close to the ones recently obtained for 
the iz-Andromedae system. A small fraction of our two-planet systems (~ 5%) ends up in the 
stability zone of the Kozai resonance. In a few cases, the triple resonance can remain stable 
for long times and a "3-D" system can form by resonant excitation of the orbital inclinations; 
such a three-planet system could be stable if enough eccentricity damping is exerted on the 
planets. Finally, in the single-planet resulting systems, which are formed when two planets 
are ejected from the system, the inclination of the planet's orbital plane with respect to the 
initial invariant plane - presumably the plane perpendicular to the star's spin axis - can be as 
large as ~ 40°. 
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1 INTRODUCTION 

The discovery of ~ 50 (June 2010) extrasolar multi-planet systems, 
whose planets have orbital characteristics quite different from the 
ones in our solar system, has opened new questions about the for- 
mation, evolution, and stability of such systems. Because of the 
lack of spatial resolution of the orbits, as well as a general belief 
that planetary systems are composed of planets on nearly co-planar 
orbits, a limited number of studies have addressed the dynamics of 
"3-D systems", namely systems composed of two or more planets, 
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whose orbital planes have a significant value of mutual inclination. 
Such "3-D" systems can be long-term stable, either following regu- 
lar secular dynamics or due to the action of some phase-protection 
mechanism, such as a mean motion resonance (MMR) or a secular 
Kozai-type resonance <Kozaill 19621 ). 

In lLibert & TsiganisI ( l2009al) . we studied the possibility that 
extrasolar two-planet systems, similar to the ones that are observed, 
can be in a stable Kozai-resonant state. Five known multi-planet 
systems that are not in MMR were selected as 'possible proto- 
types' {v Andromedae, HD 12661, HD 169830, HD 74156, HD 
155358, see Table[T]for their orbital parameters). Following a para- 
metric numerical s tudy, verified by the analytical secular theory of 
iLibert & HenrardI (l2007bl l2008l) . we found that four of these sys- 
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Table 1. Parameters of t he five exosystems an alyzed by iLibert & Tsiganisll2009ai All the parameters come from iButler et al.l ( 120061) . excepted those of the 
system HD 155358 from lCochran et al.l J2007l) . 
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terns (f Andromedae, HD 12661, HD 169830, HD 74156) are con- 
sistent with a stable Kozai-resonant state, if their (unknown) mu- 
tual inclination is ~ 45°. It should be stressed that observational 
uncertainties and/or incomplete modeling from our part (e.g. ab- 
sence of general-relativistic precession) are such that one cannot 
identify Kozai-type motion with the desired certainty. However, our 
study showed that a good fraction of the detected multi-planetary 
systems have physical/orbital characteristics compatible with the 
Kozai state. 

The results of IMc Arthur et al.l (120101) for the u-Andromedae 
system have verified that the orbital planes of planets c and d of 
this non-resonant system have a mutual inclination of ~ 30° . Thus, 
although the Kozai state is not supported by the observation, the 
large value of the mutual inclination was confirmed. Thus, the sys- 
tem is near but most likely outside the libration zone of the Kozai 
resonance. The ratio of the semi-major axes of the two planets is 
such that orbital intersections are avoided (a = ajad » 0.32) and 
so the secular dynamics are regular. 

Resonant systems of two planets can also be "3-D", since 
the orbital inclinati ons can grow under the acti on of MMRs. This 
was first shown by iThommes & Lissaueii ( 1200 3i) . who studied the 
evolution of planets trapped in a 2:1 MMR, under the effects of 
gas-driven (Type II) migrati on, inside the protoplanetary gas disc. 
iThommes & Lissaueii ( 120031) found that, once the eccentricities are 
high enough, the system that is captured in 2:1 MMR also en- 
ters an "inclination-type" resonance, which induces rapid growth 
in the inclination s of bot h planets. This work was extended by 
iLibert & TsiganisI l l2009g) . where it was shown that capture in 
higher-order resonances (such as the 5:2, 3: 1, 4: 1 and 5:1) are pos- 
sible for a wide range of migration and eccentricity damping rates 
(i.e. a and e). Moreover, it was shown that these MMRs are also 
able to excite the inclinations to high values, provided that eccen- 
tricity damping is not very strong, so that at least one of the plane- 
tary orbits has an eccentricity higher than e = 0.4. The conclusion 
of that work was that the inclination-excitation mechanism can be 
quite common in resonant systems and thus a large number of them 
may in fact represent cases of "3-D" systems. We stress here the 
fact that resonant trapping due to Type II migration is considered as 
the main formation mechanism for resonant multi-planet systems. 

Dynamical instabilities of systems with giant planets have 
been proposed to explain the orbital properties of extrasolar sys- 
tems, in particula r the highly eccentric orbits seen in many systems. 
iFord et al.l ( 120050 showed that the current orbital configuration of 



the giant planets in the u-Andromedae system probably resulted 
from planet-planet scattering with an additional planet, now lost 
from the system. A similar mechanism, accompanie d by tidal cir- 
cularization of the resulting highly eccentric orbits ( Ford & Rasid 
|2008) . is commonly invoked to explain the formation of "hot 
Jupiters" (i.e. massive planets on nearly circular orbits and with 
semi-major axes smaller than 0.03 AU). iNagasawa et al.l ( 1200 8h 
also showed that a combination of planet-planet scattering, tidal 
circularization and the Kozai mechanism in systems with three 
Jupiter-mass planets may ex plain the observed frequenci es of hot 
Jupiters (see also the work of Ipabrvckv & Tremainell2007h . 

The formation of non-resonant "3-D" systems (like the ones 
in Table [TJ - more precisely, the increase of the mutual inclina- 
tion of the orbital planes of the planets - is also thought to origi- 
nate from this violent dynamical mechanism, namely planet-planet 
scattering. This is in fact the same mechanism that is generally 
invoked to explain the eccentricity distribution of extrasolar plan- 
ets. Most of the previous studies on multi-planet scattering were 
focussed on gas-free systems composed of two planets, aiming 
at reproducing the observe d eccentri city distrib ution of extraso - 
lar planets (see for instance iFord et al. 2001. Ford & Rasioll2008l) . 
iMarzari & Weidenschillind ( l2002h explored the stability and fi- 
nal orbital properties of three-planet systems, showing that the 
most common outcome of gravitational scattering by close en- 
counters is hyperbolic ejection of one planet, the two "survivors" 
having significant val ues of eccentricity and mutu al inclination. 
IChatteriee et al.l ( |2008|) and ljuric & Tremainel (l2008h extended this 
work, showing that planet-planet scattering can reproduce quite 
well the observed eccentricity distribution. 

Moreover. lChatteriee et al.l ( l2008h have reported results on se- 
ries of planet-planet scattering simulations, where high mutual in- 
clinations (as high as 40° - 60°) between the orbits of the surviving 
planets were observed. In such simulations, the number of planets 
initially considered is typically higher than the number of surviv- 
ing planets, as the price for reaching a stable final configuration is 
typically the ejection of one (or more) planet from the system. We 
note also that Juric & Tremaine (2008) studied the systems with 
more than three planets, finding that the number of surviving plan- 
ets in their simulations was typically 2 or 3, suggesting that extraso- 
lar planetary systems are unlikely to harbor more than three giant 
planets. At present, multi-planet scattering seems to be the most 
promising scenario for forming non-resonant "3-D" systems with 
/,„„, > 30°. 
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The starting point of planet-planet scattering simulations usu- 
ally is that the systems are relatively compact when the gas dis- 
sipates, such that they can become unstable on relatively short 
time scales. Only few works investigated the combined action of 
disk torques and planet-planet sc attering (e.g. lAdams & Laughlig 
I2OO3I : iMoorhead & A dams' '2005'. for preliminary studies of two- 
planet systems). IXhommes et al. (2008) studied the formation of 
giant planets, assuming different values of the disc parameters (the 
planet-disc interactions were modeled using a A'-body code, com- 
bined with a 1-D disc model). In that work, it was shown that gas- 
driven migration tends to produce crowded systems, in which ec- 
centricity excitation due t o either resonances or pla net-planet scat- 
tering occurred. However, iMatsumura et al.l (120101) suggested that, 
according to their simulations (combining A'-body dynamics with 
hydrodynamic disc evolution), planet-planet scattering conditions 
are difficult to achieve. In fact, the existence of a large fraction 
of resonant systems suggests that unstable systems (i.e. systems in 
which planet-planet scattering sets in immediately after the gas dis- 
sipates) is not a unique result of the formation process. Of course, 
as is clear from this discussion, it is not easy to determine a priori 
what fraction of systems (and under which conditions) will reach 
a stable or a marginally unstable configuration, leading to planet- 
planet scattering. 

A possible "intermediate" evolutionary path has been noted 
recently, in studies devot ed to the origin of the solar system. 
iMorbidelli & Cridal ( l2007h studied the migration of the Jupiter- 
Saturn pair (also pairs of heavier planets) in gas discs with different 
characteristics. They showed that, at least for systems where the 
outer planet has a smaller mass than the inner one, trapping into a 
low-order MMR is accompanied by a significant drop (or even sign 
reversal) of the migration rate. Equal-mass planets (Jupiter mass) 
are also captured in resonance, but the migration direction does not 
change. Thus, the Jupiter-Satum was likely trapped in a 3:2 reso- 
nance during the gas phase and then stopped migrati ng. The results 
of this study were used by IMorbidelli et al.l ( 120071) , who demon- 
strated that the remaining outer planets of the solar system (Uranus 
and Neptune) could have followed a similar evolution, each one 
being trapped in a low-order MMR (3:2, 4:3 or 5:4) with the im- 
mediately preceding planet. Thus, at least in the case of the solar 
system, gas-driven migration can force the planets to enter into a 
multiple resonance, an analogue of the Laplace resonance in the 
Galilean moons. As noted in Morbidelli et al. (2007), a multiple 
planetary resonance is a delicate dynamical configuration; not all 
resonant ratios can be reached by all planetary masses (and mass 
ratios) and not all reso nances are long-term stable. As shown in 
IMorbidelli et alj ( 1200 7) resonant interaction can increase the plan- 
etary eccentricities, such that planet-planet scattering becomes pos- 
sible. Then, the (chaotic) dynamics will determine the orbital con- 
figuration in the final, likely "3-D", system. Moreover, even if the 
multiple r esonance remains stable, inclination excitation may oc- 
cur (as in iLibert & Tsiganisll2009bl) . again resulting into a "3-D" 
(resonant) system. 

In this paper, we wish to explore this mechanism, in the con- 
text of extrasolar systems; in particular, the formation of "3-D" 
multi-planet systems. The main point of interest in this mechanism 
is that it may apply to multi-planet systems that are not formed in 
a dynamically unstable configuration. In particular, if planets form 
on quasi-circular orbits that are sufficiently separated, the system 
can be stable on a ~ 10 My time scale, which may be long enough 
for differential migration in the gas disc to "lock" them in mul- 



tiple resonancqj. On the contrary, it assumes that migration can 
lead to the formation of a multiple resonance among the planets, 
which can (and most likely will) become unstable, due to eccen- 
tricity excitation induced by the resonance itself. Hence, the scope 
of the present paper is to examine the evolution of three-planet sys- 
tems, trapped in a triply-resonant configuration, for planet masses 
and mass ratios similar to the ones observed in non-resonant two- 
planet systems (see Table 1). We stress again that the conditions 
under which a multiply-resonant system with three (or more) giant 
planets - possibly much heavier than Jupiter - can form have to 
be investigated thoroughly in the future. The aim of our work is to 
describe the general behavior and orbital distribution of two-planet 
systems, resulting from the violent dissolution of a triply-resonant 
configuration, induced by gas-driven migration. A more detailed 
analysis, both in terms of modeling and statistics, is reserved for 
future work. Here, we present a first approach to this problem, by 
simulating the establishment of a 1:2:4 resonance and performing 
a first statistical study of the orbital characteristics of the result- 
ing systems. Let us note that the choice of this resonance seems 
in agreement with the observational data : at least the exosy stems 
HD 82943 (Beauge et al. 2008), HR8799 ( Marois et alj|2008l) . and 
Gliese 876 (.Rivera et al..2O10.) may exhibit the Laplace resonance. 
In the following section, we describe the set-up of our numerical 
experiments. The results are presented in Section 3 and our conclu- 
sions are given in Section 4. 



2 THE MODEL - NUMERICAL SET-UP 

Our initial (fictitious) systems consist of three giant planets that 
are forced to migrate such that a triple resonance can be reached. 
Then, the evolution of each system is followed, until a final, stable, 
system is formed. The masses and mass ratios of the planets were 
chosen such that they are close to the values of the non-resonant 
two-planet systems show n in Table 1. We remind the reader that 
iLibert & TsiganisI ( l2009ah found these systems to be stable for high 
values of mutual inclination. Thus, assuming the mass of the star 
to be M, = 1 Mo, we set m\ = 1.5 Jupiter masses. The mass of the 
second planet is set to mi = 1.15, 1.5 or 3 Jupiter masses, while 
m^ = mill, mi or 2 m2. Thus, in total we consider nine different 
mass configurations. Our main focus is on the semi-major axis and 
inclination distribution of those final systems that host two planets. 
The starting point in our simulations is the establishment of 
the triple resonance. Following the lines of Morbidelli et al. ( 200'5), 
we assume a two-step formation of the multiple resonance. First, 
we simulate the resonant capture of the first two planets (in or- 
der of increasing distance from the star), nji and m2. For reasons 
of simplicity we will only study capture in the 1 :2 MMR, in this 
paper The initial values of the semi-major axis of these two plan- 
ets were oi = 1 and a2 = 1.9. Our system of units is such that 
G = 1 and Mq = 1 ; thus, the period of a planet with a = I (AU) is 
T[ = In. Both initial eccentricities and inclinations are quasi-null 
(e = 0.001, i = 0.01°). All simulations are performed using the nu- 
merical integrator described in Libert & Tsiganis 2009b, in which a 
suitable Stokes-type drag force is added in the equatio ns of motion 
of the A'-body problem (foUowing lBeauge et al.ll2005r) , to simulate 
the exponential drift in semi-major axis (migration) and eccentric- 
ity (damping). The migration and eccentricity damping rates are 



' This can certainly be true for some range of masses and mass ratios. 
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Table 2. Percentage of surviving two-planet systems in our simulations of thi'ee planets in multiple resonance. For each mass configuration, the percentage of 
surviving two-planet systems with a mutual inchnation of the orbital planes higher than 40° is also shown. 
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set to l/Ta = 2.5 X 10~* and I/t^ = 5 x 10''' time units, res pec- 
tively; these values were taken from lLibert & Tsiganisll2009bl . For 
this set of parameters, the planets are captured in the resonance af- 
ter ~ 10^ Ti, when a[ x 0.9 and ai ~ 1.5. 

After mi and mi reach a stable 2:1 configuration, m^ is intro- 
duced in the simulation (with a^ = 3, 63 = 0.001 and (3 = 0.01°) 
and forced to migrate into the 1:2 MMR with mi, using the same 
recipe as above; hence the triple resonance 1:2:4 is reached. In 
this second step in each simulation, the drag terms are switched 
off in the equations of motion of mi and mi and only m^ suffers 
migration; this is done to prevent mi and mi from migrating too 
close to the star, while m^ slowly approaches the resonance. Not 
ap plying migration to m^ a nd mi also mimics the behavior shown 
in iMorbidelli & Cridall2007l where it is shown that the migration 
of two giant planets can significantly slow down once they become 
trapped in resonance, as a result of severe gas depletion in the inter- 
planetary region (formation of a "common gap"). The triple reso- 
nance is reached (but not always established, see next section) after 
another ~ 10^ Tj. The evolution of the three-planet system is then 
followed for another 5 x 10' Ti . In most cases, the system becomes 
unstable during this period, as the resonance pumps the eccentrici- 
ties of the planets. In order to improve our statistics on the orbits of 
the final systems, we increased our sample by cloning each unsta- 
ble run 10 times, taking the orbital parameters of the planets a bit 
before the instability time and adding small, random, deviations in 
the velocities of the planets. Two configurations did not go unsta- 
ble during the integration time (~ 6 x 10' Ti). For these cases, we 
cloned the final conditions of the planets and continued the integra- 
tion until the cloned systems became unstable. Thus, our sample 
consists of 90 runs in total. During the instability phase (planet- 
planet scattering) the drag terms are switched off and the evolution 
of the three-planets system is followed by a normal A'-body simu- 
lation. As expected, hyperbolic ejection of one planet is the typical 
outcome. We then focus our attention on the orbital distribution, in 
particular the mutual inclination, of the resulting two-planet sys- 
tems. 

Before we present our results some comments about our mod- 
eling procedure are in order. We fully realize that our modeling is 
far from complete. As mentioned earlier, it is not clear (and needs 
to be verified by hydrodynamical simulations) whether all possible 
three-planet configurations can lead to the formation of a triple res- 
onance. Here, we only study the 1:2:4 configuration, as the "pro- 
totype" of a triple resonance. However, the masses of the planets 
and the parameters of the disc will determine if and which triple 
resonance can be established. Concerning also the time-line of evo- 
lution, there are two more things to consider. A more gradual de- 
cay of the gas disc should probably be considered, rather than an 



abrupt switching-off of the drag terms in the equations of motion of 
the planets, during the instability phase. However, we do not con- 
sider this to be crucial, as our systems become unstable before this 
abrupt "decay" of the disc. On the other hand, no eccentricity or 
inclination damping is applied during the planet-planet scattering 
phase, which might be required for systems that reach this phase 
early enough, i.e. before the disc actually dissipates. We note, how- 
ever, that it is not easy to decide "how much" damping one should 
consider for planetary orbits with inclinations much larger than the 
width of the disc (i.e. planets moving mostly outside the gas disc); 
at least we are not aware of hydrodynamical simulations addressing 
this issue. Given the above, it is likely that our quantitative results 
are not very accurate, but we are confident that the qualitative be- 
havior is correctly reproduced. 



3 RESULTS 

Table |2] summarizes the results of our simulations. For mi = 
iMjup = m^/l, the systems became unstable before reaching a 
multi-resonant state; thus, the 1:2:4 resonance cannot be estab- 
lished for this mass configuration. This is because the planets 
are so heavy that, as m^ approached the 1:2 MMR with mi, the 
orbits become chaotic and planet-planet scattering dissolves the 
system. This was also o bserved in some of the simulations of 
iLibert & TsiganisI ( l2009bl) . In the following we will not consider 
these systems any more. 

For all the remaining mass configurations, a multiple reso- 
nance was established, according to the scheme described above. 
As the eccentricities of the planets grew, most of the systems 
evolved towards an unstable configuration, reaching eventually the 
required conditions for planet-planet scattering. The typical out- 
come was the ejection of one of the planets. The percentage of sur- 
viving two-planet systems is given in Table |2]f or each mass config- 
uration. Overall ~ 50% of all simulations resulted in a two-planet 
system. As expected, more massive planets tend to survive, while 
less massive planets tend to be ejected. Also, one of the two "sur- 
vivors" moves to a smaller orbital radius (closer to the star) while 
the other is left on a distant orbit; less massive planets are typi- 
cally scattered away from their initial orbits, while the more mas- 
sive ones remain near their initial locations. 

Hereafter we describe the orbital characteristics of the result- 
ing systems. We characterise as stable a system whose evolution re- 
mains bounded and quasi-periodic for an additional numerical inte- 
gration of 10^ Ti (or 2 X 10' Ti in case of 3-D two-planet surviving 
systems). In some cases, the secular instability of two-planet and 
three-planet systems might build up over much longer time scale. 
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Figure 1. Three-planet system in a 1:2:4 resonant configuration. The planetary masses are mi = 1.5, m2 = 1.5, 013 = 0.75M/„n. As in the lThommes & Lissauen 
l2003h mechanism, we observe an increase of the orbital inclinations for eccentricities greater than 0.2. This mechanism is described in more detail in the text. 



However, we believe that it would only reduce the number of sur- 
viving systems but not modify the general trend of the statistics 
on their orbital characteristics. In the following, we mainly focus 
on two-planet systems. We note however that one-planet as well as 
three-planet stable configurations were found. Finally, a set of 12 
runs resulted in two-planet systems with one of the planets spend- 
ing a long time on a highly-eccentric orbit (e ~ 0.9). Had we ex- 
tended the integration these systems would eventually dissolve, as 
the planets were still strongly interacting. Since we are focused on 
stable two-planet systems we decided to discard these runs. More- 
over, these runs would not provide reliable results, concerning the 
surviving one-planet systems, as our dynamical model is not com- 
plete enough to account for such long-lived, highly-eccentric orbits. 
Before presenting the results on two-planet systems a short 
note on the dynamics of 1-planet and 3-planet configurations 
should be given. 



3.1 Three-planet systems 

In the case of mild (or no) eccentricity damping, the triple resonant 
configuration typically becomes unstable, since the eccentricities 
of the planets keep increasing as migration continues. However, if 
the disc dissipates before the system dissolves, a stable configu- 
ration can be reached. Nine of our runs remain in a stable 1:2:4 
resonant configuration for the whole integration time span (and for 
an additional numerical integration of 10' Ti). What is interesting 



in these simulations is that the inclinations of the three planets in- 
crease to relatively high values. As shown in Fig.[T] the inclinations 
start growing after the eccentricities of the planets become greater 
than e ~ 0.2 and reach values of ~ 25° when one of the eccen- 
tricities exceeds ~ 0.5 - 0.6. Th is behavior is quite similar to the 
one described bv lXhommes & Lissauer (2003) , in the case of two 
planets trapped in a 1:2 MMR. Apparently, a similar mechanism is 
also active in triply-resonant configurations, such as the one studied 
here. We note that the percentage of systems following this type of 
evolution (i.e. "3-D" resonant three-planet systems) could be larger, 
had we assumed stronger eccentricity damping. 

To study the inclination excitation mechanism of three-planet 
systems in more detail. Fig. [T] describes the evolution of the crit- 
ical angles of this type of evolution. The two inner planets are 
initially captured in a 1:2 MMR, characterized by the libration of 
both resonant angles 9i = Ai - 2A2 + nri and 02 = A\ - IXj + ruj 
around 0°. It also means that the difference of the longitudes of the 
pericenter Acti2 = m^ - m2 = (Oi - 62)12 oscillates around 0°, 
i.e. the planets are in apsidal alignement (see panel e). The outer 
planet migrates until the capture in a 2/1 MMR with the second 
body at about 0.9 x 10^ yr. Only one of the two resonant angles, 
63 = A2 - 2A3 + CT2 and 64 = At - 2A3 + CT3, is in libration, which 
explains that hmi^ = mi - m^ = (O3 - 94)/2 does not oscillate in 
a first time (panel f). Thus the system is captured in a Laplace-type 
resonance, whose critical angle is cf) = Ai - 3A2 + 2A3 (panel d). 

As the three planets continue to migrate while in resonance, 
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their eccentricities increase (panel b). When their values are high 
enough, the system enters an inclination-type resonance: the angles 
Ofl = 2Ai - 4/12 + 2fii, Qp- = 2/li - 4/12 + 2^2, dp- = 2^2 - 4A^ + 2^2 

1 2 3 

and 9-2 - lAj - 4^3 + 2^3 start to librate around 180° at 4.5 x 

10' yr (panels g,h,i,j). A rapid growth of the inclinations of the 
three planets is observed, as well as the librations of the relative 
longitudes of the nodes, Af2i2 = Cli - CI2 = {Ofi -9j2)/4 and Af223 = 
^2 - fis = (^Ij^ - 6l,.2)/4 (panels k,l). 

At 6.75 X 10' yr, the system leaves the inclination-type res- 
onance, as 6,2 , 6:2 , 9-2 and 6,2 stop librating. Note that the relative 

'1 '1 '1 '1 

longitudes of the nodes fl^ and f223 still librate. The resonant an- 
gle di switches his libration center (so does Atiri2, see panel e) but 
the system remains stable due to the large values of mutual inclina- 
tions between the planets. The system still evolves in the Laplace- 
type resonance until the end of the simulation reproduced in Fig.[T] 
We expect this configuration to become unstable as the migration 
continues. Of course, the dissipation of the gas disc, sometimes af- 
ter resonance capture, could prevent the planets from this unstable 
phase and form a stable non-coplanar three-planet system. 



3.2 One-planet systems 

About 25% of our simulations (19 runs) result in one-planet sys- 
tems. Figure |2] shows such a behavior: after the ejection of a first 
planet, the remaining pair is still dynamically unstable and one of 
the remaining planets is eventually ejected from the system. In al- 
most all cases it is the heavier planet that survives. Although we 
do not have enough runs to make a proper statistical analysis, an 
interesting result is found, concerning the orbital inclination of the 
remaining planet. In 5 cases the final inclination is higher than 15° 
with respect to the initial invariant plane (the common orbital plane 
of the three planets and the presumed gas disc), which can be con- 
sidered as the plane perpendicular to the spin axis of the host star. 
The highest value of inclination recorded is 38°. 



3.3 Two-planet systems 

In ~ 50% of our runs (i.e. 40 runs), the disruption of the triple res- 
onance leads to planet-planet scattering that results in a two-planet 
system. The final systems are characterised by a large orbital sep- 
aration (a = ai/a2 small) between the two planets, large eccen- 
tricities and a large mutual inclination. An example is shown in 
Fig[3] About 30% of the systems have /„,„, > 40°, while the median 
is 30° ; this is actually the value found by iMcArthur et al.l 1 120 101) 
for the LI- And system. This result confirms our suggestion that "3- 
D" systems are naturally formed by the dynamical dissolution of 
a three-planet, multi-resonant configuration, on initially nearly co- 
planar orbits. 

The distribution of semi-major axis and eccentricity of the sur- 
viving planets can be seen in Fig. O As expected, this distribution 
is in ag reement with previous results on planet-planet scattering 
(see e.g. lChatteriee et al.ll2008r) . resembling the distribution of ob- 
served planets (although our runs are too few for a proper statis- 
tical comparison to be made). The distribution of the systems in 
terms of semi-major axis ratio, a, and mutual inclination, /„,„,, is 
given in Fig|4] As shown in this figure, 95% of the systems have 
a < 0.3. 20% of the 40 systems are located outside the 7:1, 8:1 or 
even 9: 1 MMR. The dynamical behavior of these systems is similar 
to that of the three first systems given in Table[T] i.e. v Andromedae, 
HD 12661 and HD 169830, which are outside the 5:1, 6:1 and 9:1 
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Figure 2. One-planet system obtained by the successive ejections of planet 
m2 at 71, 000 T[ and planet mi at 103, 000 T[ of the initial three-planet con- 
figuration. The planetary masses are m[ = 1.5, m2 = 1.5, m3 = 1.5-Mj,,,,. 
The inclination of the planet's orbital plane with respect to the initial in- 
variant plane - presumably the plane perpendicular to the star's spin axis - 
reaches ~ 15° in this simulation. 



MMR, respectively (we refer to lLibert & Henrardll2007al for an an- 
alytical verification of the proximity of these systems to MMRs). 
Two systems are found close to the 5:2 and 4: 1 MMR respectively, 
but no actua l resonant system is fo und. This is in agreement with 
the work of IChatteriee et al.l 1 1200 8h . where only 1% of the simu- 
lations resulted in temporary resonant capture (no dissipation was 
included in those simulations). Only two systems with high a - 
similar to the HD 155358 system of Table[T]- were found. All the 
remaining systems are practically hierarchical systems, like the HD 
74156 system in Table[T] 

Of particular interest are 1 3 of the 40 surviving two-planet sys- 
tems, which are highly non-coplanar (/,^j|, > 40°). These systems 
are in principle compatible with a Kozai resonance. Note that all 
these high-/|^„, systems were obtained when m^ was smaller than 
m\ and mi (see Table[2j. For these 13 systems, we performed an 
additional numerical integration for 2 x 10' Ti, in order to study 
their long-term stability. We found that 2 of these systems are ac- 
tually in a stable Kozai-resonant state (see Fig. [Sj, characterized 
by a coupled variation of eccentricity and inclination and libration 
of the argument of pericenter of the inner planet. Of the remaining 
systems, 6 have a stable secular behavior (not in Kozai resonance), 
while 5 systems became unstable (see Fig. |4]l. This instability is 
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Figure 3. Two-planet system formed by dynamical disraption of a triple 
resonant configuration on initially nearly co-planar orbits. The planetary 
masses are mi = 1.5, m2 = 3, 7713 = l.5Mj„p. After the hyperbolic ejection 
of the less massive outer planet, the remaining system is stable, with a large 
orbital separation between the two planets, large eccentricities and a large 
mutual inclination. 



not related to close encounters but rather to the fact that the sys- 
tems were located close to the separatrix that emanates from the 
unstable equilibrium of t he Kozai resonance, encircling the libra- 
tion zone (see Michtche nko et alj|2006l : iLibert & HenrardllJOOTbl : 
iLibert & Tsiganisii2009ai') . Let us note that Kozai resonance is not 
the only possible reason of secular instability, but in the present 
case, we have checked the five evolutions carefully and observed 
circulation of the argument of the pericenter alternately with libra- 
tion around 90° or 270° , showing the influence of the Kozai sepa- 
ratrix. 



4 CONCLUSIONS 

In this paper we examine a formation mechanism for "3-D" plane- 
tary systems, composed of two (or more) giant planets. The starting 
point of this work is the assumption that systems with three or more 
planets can be driven by Type 11 migration into a multiply-resonant 
configuration, an evolution similar to the one recently proposed for 
our solar system. This delicate dynamical configuration can then 
become unstable, not because the planet are "formed" too close to 
each other (as is typically assumed in simulations of planet-planet 
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Figure 4. Top panel: Semi-major axis (in logarithmic scale) - eccentric- 
ity distribution of the surviving two-planet systems ('square' symbols for 
inner planets, 'circle' for outer planets). This distribution seems in agree- 
ment with the extrasolar planets discovered so far Bottom panel: Semi- 
major axes ratio (in logarithmic scale) - mutual inclination distribution of 
the surviving two-planet systems. These systems typically have large values 
of semi-major axes ratio (a = aijai < 0.3), while the mutual inclination 
can be as high as 70°, with a median of ~ 30°. The 'plus' symbol denotes 
stable Kozai-resonant systems, 'cross' unstable systems and 'circle' stable 
systems on the integration time. 



scattering) but because the resonance can increase the eccentrici- 
ties of the planetary orbits, up to the point that the orbits begin to 
intersect. 

We studied a set of nine different mass configurations. We find 
that a three-planet resonance (only the 1:2:4 relation is examined 
here) can be established, as long as the outer planet is not very 
massive, in which case the system dissolves before the multiple res- 
onance is reached. In the opposite case (e.g. m^ = mill = mill), 
the resonance can even remain stable for a long time. Thus, as long 
as the disc dissipates quickly enough or exerts enough eccentric- 
ity damping on the planets as to prevent rapid eccentricity growth, 
a "3-D" multi-resonant system can form. The inclinations of all 
planets can increase to values > 20°, if one of the eccentricities 
of the planets becomes e ~ 0.5 - 0.6. As found in our runs, the 
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Figure 5. Example of a surviving two-planet system in a Kozai-resonant 
state whose dynamics is characterised by a coupled variation of the eccen- 
tricity and the inclination of the inner planet and a libration of the argument 
of the peiicenter of the same planet around 270° . 



inclinations actually start growing only when the eccentricities of 
the planets become larger th an ~ 0.2. This is reminiscent of the 
iThommes & Lissauea ( 120031) mechanism, for two planets trapped 
in a 1:2 MMR. Our results suggest that this mechanism can be also 
active on three-planet configurations. 

In ~ 90% of our runs, the triple resonance is dynamically dis- 
solved, as the planetary eccentricities grow because of the reso- 
nance, the orbits begin to intersect and the planets start encounter- 
ing each other. Form then on, the evolution is typical of what is 
observed in planet-planet scattering simulations. About 50% of the 
runs (in total) give rise to two-planet systems, while ~ 25% result in 
single-planet systems. An interesting feature of single-planet sys- 
tems is that the final orbital inclination of the planet can be as large 
as ~ 40°, with respect to the initial invariant plane (the disc plane). 

Two-planet systems, produced by this mechanism, are char- 
acterised by a large orbital separation of the two planets {a < 0.3 
in 95% of the cases); one planet is ejected on a hyperbolic orbit 
while the other two increase their orbital separation (the less mas- 
sive planet moves inwards or outwards). The outer edge of the dis- 
tribution (ttniax ~ 0.3) seems relatively sharp, although this may 
be related to poor statistics. If not, it is interesting to check in 
the future whether this outcome is related to the resonance studied 
(i.e. the 1:2:4 resonance) and whether other resonant configurations 
produce different values of o-max. About 20% of the systems have 
a a 0.3, which is near the value of planets c and d of the u-And 
system. The distribution of mutual inclinations is broad, covering 
the range - 70° . The media n is ~ 30° , which is al so the value re- 
ported for the u-And system JMcArthur et al.l2010r) . About 30% of 



the systems have /„,„! > 40°. Finally, 5% of the systems end up in- 
side the libration zone of the Kozai resonance, while several more 
systems with /„,„, > 40° become unstable not because of planet- 
planet scattering but because they are near (yet, outside) the sepa- 
ratrix that encircles the libration zone of the Kozai resonance. 

We conclude that the mechanism proposed here is quite robust 
in producing non-resonant "3-D" systems, characterised by large 
mutual inclination (median /„„„ = 30°) and large orbital separation 
(95% with a < 0.3). Hydrodynamical simulations are needed in 
order to understand the conditions under which a system of three- 
planets can be trapped in a triple resonance - and which resonant 
configurations are possible - for different mass ratios. Although we 
believe that our work captures the essential qualitative dynamics 
of the studied mechanism, more work is needed (in the sense of 
more detailed modeling and a larger sample of runs) before accu- 
rate quantitative results are obtained. We reserve this work for a 
future publication. 
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